function [BW, Y, X, Z] = detect_beads_LS_3(LS)

%这段代码的核心是利用高斯差分（DoG）来增强三维图像中的珠状结构，并结合局部最小值检测来定位这些结构。
% 它的实现包括预处理、区域掩码设定、多尺度模糊和局部特征提取等步骤，适合于检测图像中的亮度突出、呈点状的目标物。

LS = double(LS);
%figure,sliceViewer(LS);
BW = zeros(size(LS));
BW(100:end-200, 200:end-100, 19:end-19) = 1;  %指定一个感兴趣的体积范围，排除边界区域的可能噪声
BW(LS<(mean(LS(:)+100))) = 0;  %进一步排除背景噪声和低强度区域，只保留可能感兴趣的区域
LS = BW .* LS;
%figure,sliceViewer(LF);
%figure,sliceViewer(imadjustn(LF));


sigma1 = 2;
sigma2 = 4;
threshold = 0.1;
centers = gaussian_fit_3d(LS, sigma1, sigma2, threshold);

X = centers(:,1);
Y = centers(:,2);
Z = centers(:,3);


end
